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The effects of (in)homogeneity and size on the phase diagram of Lennard- Jones fluids are inves- 
tigated. It is shown that standard multifragmentation scenarios (finite equilibrated systems with 
conserved center of mass position and momentum) are implying a strong radial inhomogeneity of 
the system strongly affecting the phase diagram. The homogeneity constraint is therefore necessary 
for finite systems in order to align to the "meaning" of infinite systems phase diagrams. In this 
respect, a method which deduces the equation of state of homogeneous finite systems from the one 
corresponding to bulk matter is designed. The resultant phase diagrams show a strong dependence 
on the system's size. 
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Phase transitions have been studied from a long time 
in the limit of infinite matter. However, there are many 
physical situations away from the thermodynamic limit 
at both microscopic and macroscopic scales. Such sys- 
tems are presently of special interest as so far there is 
little knowledge about their thermodynamical behavior. 
There are thus fundamental questions related to these 
systems such as whether they do present phase transi- 
tions and, if yes, how can these be identified. (Recent at- 
tempts to extend the standard thermodynamics towards 
the "small" systems can be found in Refs. 

Highly excited nuclear systems are good "laborato- 
ries" for thermodynamics studies. Due to the van der 
Waals type of the nucleon-nucleon interaction these sys- 
tems are supposed ^ to exhibit a liquid-gas phase tran- 
sition as the classical fiuids. The connection is however 
not straight forward due to their (relatively) small num- 
ber of constituents and the presence of the (long-range) 
Coulomb force. The effect of the Coulomb interaction 
on the nuclear liquid-gas phase transitions was studied 
previously (see e.g. Refs. 0, 0|). The effects of other 
features specific to multifragmentation such as finite size 
(for attempts to address the finitness of the system within 
Hartree-Fock theories see e.g. Refs. 0,13) and degree of 
homogeneity are still mostly unknown. These aspects are 
addressed in the present paper. 

The classical Lennard- Jones 6-12 (LJ) fiuid is used for 
illustrating the paper's ideas. While neglecting for the 
quantum effects acting in nuclear matter, it has the ad- 
vantages of generality and tractability. The LJ potential 
writes: 
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Two versions of the above potential are considered herein. 
The first is the truncated and shifted (TS) LJ potential: 
v{r) = vo{r) — vo{rc) when r < r^, v{r) = when r > re- 
The second is the truncated and long range corrected 
(TLRC) LJ potential: v{r) — vo{r) when r < r^, v{r) — 
when r > Tc, corrections being subsequently included in 



order to account for the effect of the neglected tail [l^l ■ 
In the present work we use = 2.5 a. Phase diagrams of 
such fluids are subsequently constructed for various situ- 
ations. Phase transitions in finite systems being here ad- 
dressed, an adequate representaion is necessary in order 
to identify them. For example, the P(V)\t curves corre- 
sponding to an isochore canonical ensemble will prompt 
the (first order) phase transition through backbendings 
even in the case of small systems 0, |^. Pressure can 
be easily evaluated starting from its canonical definition: 
P = TdliiZij3,V)/dV, where Z{f3,V) is the system's 
canonical partition function: 
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where (3 =1/T and At is the thermal wavelength. The 
result is the following virial expression: 
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where ( ) has the meaning of canonical average. When 
conservation of the system's center of mass (cm.) and 
cm. momentum is also considered in Z{l3, V), one simply 
has to change A into A — 1 in eqs. 10) [the power of At] 
and (jSJ [the first term] . While rather redundant for very 
large systems {A oo), these constraints give important 
effects when dealing with small systems. Expressions like 
© can be easily evaluated by means of Metropolis sim- 
ulations. Knowing the formal expression of Z{P, V), 
the statistical simulation of the corresponding canonical 
ensemble at constant volume is straight-forward: Pick a 
randomly chosen particle and make a random variation 
of its initial position (in a fixed interval AV^) . Then, con- 
sider this move according to the acceptance exp{—PAvt), 
where Avt is the change in the system's total potential 
energy. The simulation can be easily adapted to the case 
of the conserved system's cm. and cm. momentum: 
Two randomly chosen particles are moved simultaneously 
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with Ar and — Ar such that the system's cm. remains 
fixed (Ar being a variation randomly chosen in the vol- 
ume interval AV^) . Then, the move is considered accord- 
ing to the same acceptance. 

The present study is started by considering the case of 
a system at fixed temperature T, composed of A particles 
interacting via TS LJ potential, contained into a spher- 
ical recipient of volume V, having the cm. constrained 
to coincide with the center of the recipient. The cor- 
responding canonical ensemble is simulated as described 
earlier. Then, isothermal P{V) curves can be evaluated 
by means of eq. lj2Jl, adapted to the considered conser- 
vation laws. Subsequently, the borders of the liquid-gas 
coexistence region are evaluated by performing Maxwell 
constructions on all P{V)\t curves bending backwards. 
The resulting phase diagrams (in temperature vs. density 
representation) corresponding to two systems of different 
sizes A ~ 20 and yl = 50 are represented in the upper 
part of Fig. 1. The increase of the critical temperature 
with the system's size can be observed. A striking fea- 
ture of these phase diagrams is the small densities cor- 
responding to the borders of the liquid-gas coexistence 
regions. In this respect, note that the liquid border is 
situated at densities smaller than 0.2 ct"'^ which differs 
a lot from the LJ phase diagrams corresponding to infi- 
nite homogeneous systems jinl | where the liquid border 
goes up 0.8 a~^. The reason for this discrepancy can be 
easily understood from the lower part of Fig. 1. There, 
the radial density profiles corresponding to three sam- 
ple points chosen from the coexistence line of the A=50 
system are represented. In all three cases the system 
appears to be inhomogeneous, its density varying from 
higher values (towards the center of the recipient) to very 
small ones (towards the recipient's walls). The low den- 
sity tails of these radial profiles are therefore inducing 
the above mentioned effect on the global density of the 
system (calculated as A/V). In particular, note that at 
small values of the distance from the recipient's center 
the densities of the considered sam ple p oints tend to be 
consistent with the ones from Ref. [lO|. Obviously, this 
inhomogeneity effect is dictated by the cm. conserva- 
tion constraint which forces the larger clusters to stay 
towards the center of the recipient and the smaller ones 
towards the borders. This example is particularly im- 
portant in the context of multifragmentation where the 
system's cm. is naturally conserved from event to event 
and similar density profiles are expected to occur. 

Therefore, in order to meaningfully compare the phase 
diagram calculated for finite systems with the ones cor- 
responding to infinite systems an extra constraint has 
to be imposed: the homogeneity. The infinite and ho- 
mogeneous system can be simply approached within the 
(above-described) Metropolis simulation by implement- 
ing a cubic recipient with periodic boundary conditions 
(PBC). (No cm. conservation constraint has to be im- 
posed.) A number of 200 particles interacting via TLRC 
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FIG. 1; Upper panel: Phase diagrams for two systems [A — 
20 and A = 50) with conserved cm. position and cm. mo- 
mentum. Three sample points from the A = 50 phase diagram 
are represented by squares. Lower panel: Radial density pro- 
files corresponding to the sample points from the upper part 
of the figure. The correspondence between curves and sample 
points is given by numbers. 



or TS LJ potentials are placed in a cubic box with PBC. 
The very good agreement between our calculated phase 
diagram corresponding to the TLRC case (see Fig. 3, the 
"oo" curve) and the one from Ref. is confirming the 
accuracy of the method here employed. The PBC sim- 
ulations are used for evaluating the bulk "virial energy" 
per particle, defined as: 
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Eq. © can be now translated into: 
P = p{T-V), 



(4) 
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where p is the system's density. Therefore, in an infinite 
and homogeneous system, the virial energy per unity of 
volume can be written as: V = pV. Given the above 
definitions one may perform surface corrections for eval- 
uating the V term corresponding to finite systems. Since 
the system is supposed to be homogeneous, its total virial 
energy can be expressed as the difference between a bulk 
term and a surface one: 



Vt - Vb{A - AJ) ^ %{V - Vsf), 
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where the index b specifies that the respective term is a 
bulk one, As and Vg are respectively the number of par- 
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FIG. 2: Pressure versus density curves corresponding to sys- 
tems of various sizes (^4) at fixed temperature, T — 1.05e, 
evaluated with eq. JSJ . The corresponding size of the system 
is specified on each curve. 

tides from the surface and the "volume" of the surface, 
and finally, 1 — / is the ratio between the virial energy 
of a particle from the surface and the virial energy of a 
bulk particle. While for very large values of the recipi- 
ent's radius {R) 1 — / is rigorously equal to 1/2 (i.e. the 
particles on the surface have half the number of nearest 
neighbors they have in bulk), for small values of R (the 
case of small systems) curvature corrections have to be 
applied to this factor. In the spirit of the previous defi- 
nition, / can be fairly approached by the ratio between 
the surface of spherical cap situated outside the recipient, 
corresponding to a sphere of radius a having the center 
on the surface of the recipient and Aira'^ . After some 
algebraic manipulation one gets: 

Note that when R ^ oo then / ^ 1/2. Considering the 
surface "width" equal to a (i.e. no two particles from 
that region, situated on the same radius, attract each 
other) one can re-express eq. © as: 




Eq. lO can be further expressed in terms of p using the 
identity a/R= [iir p /(3A)]i/3. Subsequently, eq. © 
[with V given by eq. ©] can be applied for calculating 
the pressure of a finite system of size A at various values 
of p. 
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FIG. 3: Phase diagrams corresponding to the TLRC LJ fluid, 
corresponding to different sizes of the system. The system's 
size is specified on top of each diagram. Points calculated 
via Maxwell constructions are represented with symbols. Full 
lines are fits of eq. IIUII on the calculated points. 

An illustration of the results of the method is given in 
Fig. 2 where P{p)\t curves, corresponding to the tem- 
perature r = 1.05 e are represented for various sizes of 
the system. One can observe that the depth of the back- 
bending of the curve is diminishing as the size of the 
source is decreasing such that at system sizes as small as 
A = 20 the backbending completely disappears. 

Maxwell constructions have been further performed on 
the P{V)\t curves corresponding to various values of T 
and A allowing thus the construction of phase diagrams 
for various sizes of the system. The results corresponding 
to the TLRC and TS LJ forms of potential are given in 
Figs. 3 and 4 respectively. The points from the coexis- 
tence region borders obtained via Maxwell constructions 
can be further interpolated in order to get an estimation 
of the critical point and a clearer view on the systems' 
phase diagrams by means of the Guggenheim scaling re- 
lations for the coexistence curve: 

p±= Pc{l + ae±htP), (10) 

where p+ corresponds to the liquid branch, p_ to the va- 
por branch of the coexistence curve, pc is the critical den- 
sity, e = (Tc — T)/Tc {Tc being the critical temperature) 
and (3 is the critical exponent of the coexistence curve. 
Eq. I|10|) is fitted on the points obtained via the Maxwell 
construction method by adjusting the parameters: Tc, 
Pc /3, a and b. As observed in Figs. 3 and 4 the quality 
of the obtained fits is very good. (It is worth mentioning 
that for the TS potential, infinite system case, j3 = 0.34 
which sharply corresponds to the liquid-gas universality 
class.) For both considered potentials the critical tem- 
perature is drastically decreasing and the critical density 
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FIG. 4: Same as Fig. 3 but for a TS LJ fluid. 



This method makes possible the deduction of the P{V)\t 
curves corresponding to any size of the system, A, by us- 
ing the information embedded in the bulk P{V)\t curves. 
The resulting phase diagrams corresponding to homoge- 
neous systems of various sizes show a strong dependence 
on the size of the system. In this respect, the critical tem- 
perature is drastically decreasing when the size of the 
system is reduced. For example, important deviations 
from the bulk phase diagram can be noticed even for 
systems as large as A — 10^. This means that the largest 
equilibrated systems formed in nuclear multifragmenta- 
tion experiments are not even close to the bulk limit. 
To this effect one should, of course, couple the inhomo- 
geneity (discussed earlier) and the Coulomb ones. And 
as shown in Ref. Q Coulomb is independently bring- 
ing a very important contribution towards lowering the 
system's critical point. 

This work was supported by the Alexander von Hum- 
boldt Foundation. 



is increasing with decreasing the size of the source. In 
particular, note that important deviations from the "oo" 
phase diagram can be observed even for the phase dia- 
gram of a system as large a.s A = 10^ (see Fig. 4). This 
result is particularly important since it proves that ther- 
modynamically the bulk limit is not (even by far) attain- 
able in nuclear multifragmcntation experiments (where 
the equilibrated systems formed are usually smaller then 
A = 300). 

Summarizing, homogeneity and size effects on the sys- 
tem's phase diagram have been discussed in the frame- 
work of the classical LJ fluid. Phase diagrams are con- 
structed using a method adequate for revealing the co- 
existence region even for small systems: Maxwell con- 
structions on canonical P(V)\t curves. It was shown 
that standard scenarios like finite systems with con- 
served cm., physically consistent with the multifragmen- 
tation phenomenon imply a strong radial inhomogene- 
ity in the systems importantly affecting their phase di- 
agrams. Thus, in order to align to the meaning of the 
phase diagrams corresponding to the infinite fluids 
an extra constraint has to be introduced for the case of 
the finite systems: homogeneity. To this aim, a method 
for constructing the P{V)\t curves corresponding to ho- 
mogeneous finite systems based on making surface cor- 
rections on the bulk virial energy term was designed. 
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